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Abstract. We report on the results of the spectral and timing analysis of a BeppoSAX observation of the microquasar system 
LS 5039/RX J1826.2-1450. The source was found in a low-flux state with Fx(l-10 keV)= 4.7 x 10"'^ erg cm"- s"', which 
represents almost one order of magnitude lower than a previous RXTE observation 2.5 years before. The 0. 1-10 keV spectrum 
is described by an absorbed power-law continuum with photon-number spectral index T = 1.8 ± 0.2 and hydrogen column 
density of A'h = l-O^JJ* x 10^^ cm"^. According to the orbital parameters of the system the BeppoSAX observation covers the 
time of an X-ray eclipse should one occur. However, the 1.6-10 keV light curve does not show evidence for such an event, 
which allows us to give an upper limit to the inclination of the system. The low X-ray flux detected during this observation is 
interpreted as a decrease in the mass accretion rate onto the compact object due to a decrease in the mass-loss rate from the 
primary. 
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1. Introduction 

LS 5039/RX J 1826.2- 1450 was first identified as a new mas- 
sive X-ray binary by Motch et al. (11997*. Interest in this source 
has grown significantly because it has turned out to be a source 
of persistent radio-emitting relativistic jets, and because it is 
believed to be physically associated with the y-ray source 
3EG J1824-1514(Paredes et al 2000l l2002t . 

The first radio detection was reported by Marti et al. ( 1998 1 
using the Very Large Array, but the discovery of radio jets was 
only possible when the source was observed at milliarcsecond 
scales with the Very Long Baseline Array (Paredes et al. 2000 1. 
The radio emission is persistent, non-thermal and variable, but 
no strong radio outbursts or periodic variability have been de- 
tected (Ribo et al. [19991 hereafter R99; Rib6 l2U02t . 

In the optical band LS 5039/RX J 1826.2- 1450 appears as 
a bright y=11.2, 06.5V((f)) star showing little variability on 
timescales of months to years (Clark et al. 1200 U . Variations 
of ~0.4 mag have been reported in the infrared (H and K 
bands) but no obvious mechanisms for such variability have 
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been proposed (Clark et al. "2001'!. Recent studies show that 
LS 5039 is a runaway system moving away from the Galactic 
plane with a total systemic velocity of ~150 km s"' and a 
component perpendicular to the Galactic plane larger than 
100 km s-i (Ribo et al. 2002 , McSwain & Gies 2002 ). The or- 
bit ofLS 5039/RX J 1826.2- 1450 was first studied by McSwain 
et al. (2001 1, who found that LS 5039/RX J1826.2-1450 is a 
short-period, Porb = 4.117 + 0.011 day and highly eccentric, 
e = 0.41 ± 0.05, system. Spectroscopic observations carried 
out in July 2002 by Casares et al. ('2003^ seem to confirm the 
orbital period of the system, but neither confirm nor reject the 
proposed eccentricity. 

In the X-rays LS 5039/RX J1826.2-1450has been studied 
by R99 using RXTE data. The X-ray timing analysis indicates 
the absence of pulsed or periodic emission on time scales of 
0.02-2000 seconds. The source spectrum is well represented 
by a power-law model, plus a Gaussian component describ- 
ing a strong iron line at 6.6 keV. Significant emission is seen 
up to 30 keV, and no exponential cut-off at high energy is re- 
quired. The X-ray luminosity has been shown to be compatible 
with accretion from the stellar wind of the optical companion 
(McSwain & Gies UM^ . 

One of the open questions that remains to be solved is 
the nature of the compact object of the system as no conclu- 
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Fig. 1. MECS background subtracted light curves of 
RX J 1826.2- 1450 in different energy ranges. Each point rep- 
resents 60 minutes. Time zero corresponds to JD 2,451,825.75. 
The expected time of periastron (P) and inferior conjunction 
(IC) of the primary are indicated in the top panel, where the 
error bars represent the uncertainty in the zero point of the 
ephemeris. The time difference between these two points 
is 2.2 X 10"* seconds. An estimate of the unabsorbed X-ray 
luminosity according to the model of Table ^ is also given in 
the right axis of the top panel. 



sive data exists. The absence of X-ray pulsations and exponen- 
tial cut-off at ~30 keV and the persistent non-thermal emis- 
sion favours the black-hole classification (R99). On the other 
hand, the value of the mass function of the system /(m) = 
0.00103 M0 (McSwain et al. 1200 11 1 and the characteristics of 
its accretion mechanism (McSwain & Gies 2002jl appear to be 
more consistent with a neutron star companion. 

In this work we present an analysis of the X-ray timing and 
spectroscopic properties of RX J 1826.2- 1450 using the nar- 
row field instruments on board the BeppoSAX satellite. This 
observation permits us to derive a more accurate value of the 
absorption to the source and study its long-term X-ray variabil- 
ity. We also present optical spectroscopic observations around 
the Ha line carried out at the Skinakas Observatory. 

2. Observations and results 

2A. X-ray observations 

LS 5039/RX J1826.2-1450 was observed by BeppoSAX on 
October 8, 2000 (JD 245 1825.79-245 1826.71) for about 80 ks. 
BeppoSAX (Boella et al. 19_97a ) carried four narrow field in- 
sti-uments (LEGS, MECS, HPGSPC and PDS) and a wide 
field camera (WFC). Here we report on the results from the 
lower energy detectors: the Low-Energy (LEGS) and Medium- 
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Fig. 2. LEGS energy spectrum of RX J 1826.2- 1450 and the 
best-fit model displayed in Table[3 



Energy (MEGS) Concentrator Spectrometers. They consisted 
of a set of four and three X-ray concentrators respectively to- 
gether with imaging gas scintillation proportional counter de- 
tectors located at the focal planes. In the LEGS (Parmar et al. 
I1997L one of these detectors was sensitive to X-rays in the en- 
ergy range 0.1-10 keV, while the other three covered an energy 
range of 1 .3-10 keV. The MEGS (Boella et al. ll997b> operated 
in the energy range 1.6-10 keV. Each detector unit provided a 
total collecting (geometrical) area of ~ 124 cm^, an energy res- 
olution of < 8.8 % at 6 keV and a maximum time resolution of 
16yL(s. Note that one of the MEGS became inoperative on May 
9, 1997 due to a fault in the unit's gas cell high voltage supply. 

Light curves and spectra were extracted from circular re- 
gions with a radius of 4' for both the LEGS and the MEGS 
instruments. We selected two source-free regions (of radius 4' 
each) in the field of view, and obtained a mean light curve in 
order to account for the instrumental background. For the spec- 
tral analysis, the background was obtained from a blank field 
using a region similar in size and position to the source extrac- 
tion region. 

The X-ray intensity in the energy range 1.6-10 keV de- 
creased with time from a value of ~0.06 MEGS count s"' at the 
beginning of the observation to ~0.03 MEGS count s ' at the 
end (Fig.0. A fit to a constant gave an unacceptable reduced 

(>2). However, adding a linear term with negative slope we 
obtained a reduced ~\. An F-test gave a probability of less 
than 0.7% that this improvement in the quality of the fit could 
occur by chance. The flux variability occurs at higher energies. 
Below ~3 keV the flux remains fairly constant (Fig.[0. 

Acceptable fits of the X-ray energy spectra were obtained 
with an absorbed power-law model for both LEGS (Fig.|2} and 
MEGS (Fig. Table[2shows the best-fit spectral parameters. 
The power law was modified by a low-energy absorption com- 
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Fig. 3. MECS energy spectrum of RX J 1826.2- 1450 and the 
best-fit model displayed in Table [0 

ponent of the form e-^no-(£)^ where A^h is the equivalent hydro- 
gen column density and cr{E) the photoelectric cross-section 
of Morrison & McCammon ( 1983 1. The power-law index and 
the hydrogen column density have values of ~1.8 and ~ 1.0 x 
10^^ cm"^, respectively (see discussion in Sect. 13. 1> . The LECS 
(0.1-5 keV) and MECS (1.6-lOkeV) absorbedX-ray luminos- 
ity was 1.9 X 10^^ erg s"' and 4.6 x lO^"' erg s"', respectively, 
assuming a distance of 2.9 kpc (Ribo et al. l2002t . 

2.2. Optical observations 

Optical spectroscopic observations were made with the 1.3 m 
telescope of the Skinakas Observatory in Crete (Greece) on the 
night of July 26, 1999 and September 10, 2002. The instru- 
mental set-up consisted of a 2000x800 ISA SlTe CCD and a 
1302 1 mm"' grating, giving a dispersion of ~1 A/pixel. The re- 
duction of the spectra was made using the STARLINK Figaro 
package (Shortridge et al. 2001 1, while their analysis was per- 
formed using the STARLINK Dipso package (Howarth et al. 

The wavelength coverage of the observations included the 
Ho- line (/i6563 A). The strength of the stellar wind emis- 
sion in Wa provides information about the mass-loss rate in 
O-type stars (Puis et al. 1996 and references therein). The mea- 
surement of the equivalent width of the Ha line (EW(Ha)) in 
LS 5039 is hampered by the presence of other stellar Unes, es- 
pecially He II /16527. Nevertheless, when the emission is low 
the main factor contributing to the uncertainty in the equivalent 
width is the difficulty in obtaining the intensity of the contin- 
uum. In Table 12 we give the values of EW(Hq') as contempora- 
neous as possible with the X-ray measurements. 



3. Discussion 

3.1. Hydrogen column density 

The RXTE observation reported by R99 was not able to pro- 
vide a reliable value of the absorption to the source due to 
the insensitivity of the detectors to X-rays below 2.5 keV, 
where most of the interstellar absorption of X-ray photons takes 
place. The sensitivity of the BeppoSAX narrow field instru- 
ments to X-rays below 2 keV, especially LECS, allows us to 
derive a more accurate value of the hydrogen column density, 
A^H = ^-^-ot ^ 10^^ cm-2. Assuming E(B - V) = Av/3.1 
(Rieke & Lebofskv I1985> the X-ray absorption can be con- 
verted into optical extinction by means of the empirical rela- 
tions of Predehl & Schmitt (.1995j . to obtain E(B-V) = 1.8+°-^, 
or those of Ryter et al. J1975> and Gorenstein J1975t . to obtain 
E{B -V)- 1.5^04. These values of E{B - V) are considerably 
higher than the unreliable value of E(B -V) * 0.4 derived from 
the PCA/RXTE observation (R99). 

On the other hand, optical photometry of LS 5039 can be 
found in Drilling ( 1991 1 and Lahulla & Hilton (il992j who gave 
a value for the colour (B - V) - 0.95 + 0.01 and Clark et al. 
12001 ), who reported (B-V)^ 0.85 + 0.02. Using the aver- 
age of these values and assuming an intrinsic colour index of 
(B - V)o = -0.30 ± 0.02 for an 06.5V star (Lejeune & Scharer 
I20(yn we obtain a colour excess of E{B - V) = 1.2 ± 0.1. The 
extinction obtained from the BeppoSAX observation is some- 
how higher than (although statistically consistent with) that de- 
rived from the optical observations. One explanation for this 
diff'erence might be the fact that the hydrogen column density 
is more sensitive to absorption by both interstellar and intrin- 
sic (i.e. in the vicinity of the X-ray source) matter, while the 
optical indicators would mainly account for the interstellar ab- 
sorption. The presence of significant amounts of cold matter in 
the surroundings of the system is supported by the detection of 
the relatively strong iron (equivalent width of 0.75 ±0.06 keV) 
emission line at 6.6 keV in the 3-30 keV X-ray spectrum of 
LS 5039/RX J 1826.2- 1450 (R99). Unfortunately, due to lim- 
ited signal-to-noise ratio no line is detected in the BeppoSAX 
data. 

3.2. Absence of X-ray eclipse 

Using the ephemeris of McSwain et al. (I2001t . the BeppoSAX 
observation covers the orbital phase from 0.89 to 0.11. Note 
that this observation took place just four days after the pe- 
riastron date used by these authors in their ephemeris, i.e., 
JD 2,451,822. 12+0.09. Consequently, our determination of the 
orbital phase is not strongly affected by the lack of accuracy 
in the value of the orbital period (Porb = 4.1 17 ± 0.01 1 days). 
Hence, periastron took place on JD 2,45 1,826. 24±0. 09, or at 
4.2 ± 0.8 X 10"* s in Fig. ^ The epoch of inferior conjunc- 
tion of the optical star, which would correspond to the time 
of an X-ray eclipse, is 0.25 days after periastron or at or- 
bital phase 0.06 (McSwain et al. 2001 1, which corresponds to 
JD 2,45 1,826.49+0.09, or to 6.4 + 0.8 x lO" s in Fig. [T] The 
duration of a possible eclipse would last less than 0.2 days 
(1.7 X 10"* s) for Mont = 40 Mq, /^opt = 10 Rq and ; = 90° 
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Table 1. Spectral fit results for the absorbed power-law model. 
Uncertainties are given 68% confidence for one parameter of 
interest. 



Table 2. X-ray flux and Ha equivalent width measurements of 
LS 5039. Errors in EW(Ha) are <10%. 



Parameters 




Best-fit values 


LEGS 






1.0!«txl022 


r 




1.8±6.3 


Fx(0.1-5keV) (erg 


cm ^ s ') 


1.9 X 10^'- 


(dof) 




1.25 (12) 


MECS 


A^H (cm-2) 




1.2±0.4xl02- 


r 




1.8±0.1 


Fx (1.6-10 keV) (er 


J cm ^ s ' ) 


4.6 X 10-'2 


X'r (dof) 




0.87 (17) 



(shorter for a lower inclination and an eventual lower mass 
for the primary) and hence it should have been detected by 
BeppoSAX. 

Although the light curve does not show statistical signifi- 
cant evidence for an eclipse, the variability increases between 
periastron and inferior conjunction, possibly indicating some 
kind of distortion of the stellar wind when the compact star 
finds itself at minimum distance from the primary. There is also 
a trend for the count rate to decrease as the source approaches 
inferior conjunction. This trend is detected at energies above 
3 keV, but not at lower energies, contrary to what would be ex- 
pected if this was due to a change in absorption, as it happens 
in Cygnus X-1 (Wen et al. I1999> . Therefore, this decrease in 
the hard X-ray flux during the BeppoSAX observation is prob- 
ably intrinsic, and may reflect changes in the accretion rate due 
to orbital motion. 

The lack of eclipses allows us to set an upper limit on the in- 
clination of the system /. This impUes, in the case of a spherical 
companion (no photometric ellipsoidal modulation has been 
found to date), that cos / > Ropt/r, where /?opt is the radius of the 
optical companion an r is its distance to the compact object at 
inferior conjunction (phase 0.06). Assuming a mass of 40 
and a radius of 10 R© for an 06.5V star (Howarth & Prinja 
I1989> . the canonical mass for a neutron star (1.4 Mq) and the 
orbital parameters of McSwain et al. ( 2001 1, e = 0.41 ±0.05 and 
^'oib =4.1 17±0.01 1 days, we obtain an upper limit for the incli- 
nation of the system of / < 66 ± 2°, where the upper and lower 
values correspond to the range of eccentricities given above. 
This incHnation is consistent with the one, //, derived from the 
mass function, f(m) = 0.00103 ± 0.00020 Mq (McSwain et al. 
1200 1>. for the same parameters, namely, // = 60°. 

Note that OB stars in X-ray binary systems may be less 
massive than the corresponding field stars of the same spectral 
type (Kaper ^OOIV A decrease of the mass of the optical com- 
panion by about 10% reduces // to 54°, while the upper limit 
hardly changes. On the other hand, assuming Mopt = 40 M© 



Mission 



Date 

X-ray obs. 



EW(lla) Date 
(A) Opt. obs. 



RXTE 8/02/1998 40* 

ASCA 4/10/1999 10.9* 

BeppoSAX 8/10/2000 4.9* 

CHANDRA 10/09/2002 10.0'' 



2.2'-- 
2.8* 
3.1'' 
2.8* 



Aug 1998 
26/7/1999 
Oct 2000 
10/9/2002 



xlO erg cm ^ s 



and the above mentioned orbital parameters Ribo et al. J2002t 
derived an upper limit for the mass of the compact object of 



' in the energy range 0.3-10 keV. 
* This work. 
McSwain & Gies <2002l . 
Miller (private communication). 



9 Mq (in order to avoid breackup speed of the primary once 
Vrot sin / is known). For a 9 M© black-hole companion a similar 
analysis yields ; < 68 ± 2° and // = 8.5°. Similar values are 
obtained if we decrease the mass of the optical companion by 
about 10%. In either case, the lack of X-ray eclipses constrains 
the inclination of the system to values smaller than 70°, which 
is clearly compatible with those found using the mass function 
of the system. 



3.3. Long-term X-ray variability 

In order to study the long-term X-ray variability of 
LS 5039/RX J 1826.2- 1450 we obtained the X-ray luminos- 
ity in the same energy ranges as given in previous missions. 
The value of the ROSAT 1996 October 0.1-2.4 keV (Motch 
et al. 19971 and the RXTE 1998 February 3-30 keV (R99) X- 
ray luminosities are 7.1 x lO-'-' erg s"' and 5.2 x lO-'"^ erg s"', 
respectively. The 0.1-2.4 keV X-ray luminosity obtained with 
LEGS is 6.6 x 10^^ erg s ' and that obtained with MECS in the 
energy range 3-30 keV is 8.0 X 10^^ erg s"' (derived assum- 
ing the model given in Table This observation took place 
in 2000 October All values are for an assumed distance of 
2.9 kpc (Ribo et al. 2 002. ). Clearly the source was in a low-flux 
state during the BeppoSAX observation (about one order of 
magnitude fainter) when compared to the epoch of the ROSAT 
and RXTE observations. As explained above this low flux, 
however, cannot be accounted for by invoking a partial X-ray 
eclipse. We have also inspected archival X-ray data and found 
an ASCA observation of LS 5039/RX J 1826.2- 1450 made on 
1999 October 4. We have obtained the X-ray flux in the en- 
ergy range 0.3-10 keV, shown in Table|2] which is between the 
RXTE and BeppoSAX measurements. In addition, a similar X- 
ray flux has been recendy detected with CHANDRA (Miller, 
private communication), and is also reported in Table|2 

We also studied the long-term optical variability of the 
primary by measuring the strength of the Ha line. Table |2l 
shows the values of the Ha equivalent width together 
with the X-ray flux from the missions that have observed 
LS 5039/RX J1826.2-1450. The Ha emission was weak- 
est during the BeppoSAX observation in October 2000 and 
strongest a few months after the RXTE pointing in February 
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1998, when the X-ray emission was also strong (since the line 
is in absorption, the smaller EW(lia) the stronger the emis- 
sion). Moreover, very similar EW{Ha) are found contempora- 
neous to the ASCA and CHANDRA observations, as the de- 
tected X-ray fluxes are, and also between the EW(lia) mea- 
surements contemporaneous to the RXTE and BeppoSAX ob- 
servations. Hence, there is a trend in the sense that the lower the 
emission in Ha, the lower the observed X-ray flux. All these re- 
sults suggest that a decrease in the mass-loss rate from the pri- 
mary induces a decrease in the emitted X-ray flux of the binary 
system. In this context, we propose that changes in the mass- 
loss rate from the primary, reflected in the change of EW{lia), 
can account for a change in the mass accretion rate onto the 
compact object, and hence explain the detected changes in X- 
ray luminosity. 

On the other hand, since this binary system has a high ec- 
centricity, we cannot rule out a priori that the detected vari- 
ations in the X-ray flux are not due to orbital variability. If 
the mass transfer in LS 5039/RX J 1826.2- 1450 occurs via the 
strong stellar wind of the optical primary then, given the high 
eccentricity of the system, some variations in the X-ray lumi- 
nosity with the orbital phase should be expected. The fraction 
of the stellar wind captured by the compact object depends on 
its velocity relative to the wind. Also, variations in the den- 
sity of the wind may produce changes in the amount of ab- 
sorption material throughout the orbit giving rise to the X-ray 
flux variations. Orbital flux variations have been seen in many 
high-mass X-ray binaries powered by stellar wind accretion: 
4U 1700-37 (Haberl & Dav [T992t . GX 301-2 (Haberl [T99Tt . 
Vela X-1 (Haberl & White mU), 4U 1907-1-09 (Marshall & 
Ricketts .l980j . With its small orbit and early-type optical pri- 
mary LS 5039/RX J1826.2-1450 is expected to share many 
similarities with these systems. 

However, a purely Bondi-Hoyle accretion model (Bondi 
I1952> does not reproduce the observations. Figure 0] 
shows the expected variability of the X-ray luminosity in 
LS 5039/RX J 1826.2- 1450 using a spherically symmetrical 
Bondi-Hoyle accretion model. For each orbital phase and using 
the formulae given in the Appendix, we computed the binary 
separation, the orbital and wind velocities, the relative velocity, 
the wind density, the accretion rate and finally the X-ray lumi- 
nosity. We note that this is a bolometric unabsorbed X-ray lu- 
minosity, to be compared to the much lower observed luminos- 
ity in the range 1.6-10 keV (top panel of Fig.Q. Although this 
discrepancy could also be caused by the lack of precise knowl- 
edge of the parameter values used in the model, there are two 
strong arguments against this model. First, we should expect a 
variation in flux of at least a factor of 2 within the phase inter- 
val of the BeppoSAX observations as follows: starting with a 
gradual rise, peaking smoothly around phase 0.02 (~ 5 x lO'* s 
in Fig.[0 and ending with a gradual decrease. Second, since the 
BeppoSAX observation took place during periastron the X-ray 
luminosity should be the largest, contrary to what it is observed 
(see Table|5J. 

Although it is clear that detailed observations covering a 
full orbital period are needed in order to better study and model 
the orbital X-ray variability of LS 5039, it seems unlikely that it 
can be responsible for the variations of one order of magnitude 



2.0x10' 




-0.5 0.0 0.5 1.0 1.5 
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Fig. 4. Expected variability of the X-ray luminosity in LS 5039 
due to the significantly eccentric orbit of the compact compan- 
ion around the 06.5V((f)) star with a period Porb - 4.117 d, 
using a spherically symmetrical Bondi-Hoyle accretion model 
and the parameters quoted in the text. The thick line corre- 
sponds to e = 0.41, the nominal eccentricity value in the 
McSwain et al. ^2001* solution, while the lower and upper thin 
lines are for e - 0.36 and e - 0.46, which correspond to the 
error limits quoted by the authors. Phase corresponds to peri- 
astron, while the shaded areas correspond to the phase interval 
covered by our BeppoSAX observation. 

quoted in Table |2l Therefore, we conclude that this variability 
may be due to changes in the stellar wind mass-loss rate on 
longer timescales. 

4. Conclusions 

The main conclusions of the analysis of our BeppoSAX obser- 
vation of the microquasar LS 5039/RX J 1826.2- 1450 can be 
summarized as follows: 

1 . We have obtained the most accurate value of the hydrogen 
column density to date: A^h = l-0^o 3 x 10^^ cm"^. The de- 
rived colour excess is compatible but higher than the one 
obtained from optical photometry alone, suggesting that 
there might be intrinsic absorption within the binary sys- 
tem orbit. 

2. Our observation took place near inferior conjunction, and 
no X-ray eclipse has been detected. This fact has allowed 
us to compute upper limits for the inclination of the binary 
system for different masses of the compact object: / < 66° 
for Mx = 1.4 Mq and / < 68° for Mx = 9 Mq. These val- 
ues are compatible with the ones obtained using the mass 
function. 

3. We have found hints of a correlation between the mea- 
sured EW(lia) and the observed X-ray flux on timescales 
of years, in the sense that the stronger the optical emission 
the higher the X-ray flux. In this context, we have suggested 
that a decrease in the mass-loss rate from the primary in- 
duces a decrease in the mass accretion rate onto the com- 
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pact object, and hence a decrease in the observed X-ray 
flux. 

4. A simple Bondi-Hoyle accretion model (without taking 
into account a varying absorption) cannot explain the X- 
ray variability because it predicts an increase in the X-ray 
flux near periastron that is not observed. This may happen 
because the model is very sensitive to some parameters not 
yet very well constrained, like the eccentricity and the wind 
velocity. However, the fact that no significant variability is 
detected near periastron suggests that the observed variabil- 
ity in the X-ray flux on timescales of years is due to changes 
in the mass-loss rate of the primary. An X-ray monitoring 
during a full ~4 d orbital cycle would clarify these ideas. 

Appendix: Orbital Bondi-IHoyle X-ray variability 

An order of magnitude estimate of the orbital X-ray variability 
can be obtained using the Bondi-Hoyle accretion model. The 
accretion luminosity is then computed by means of the simple 
formula: 



GMxM 
Rx ' 



(1) 



where G is the gravitational constant, Mx and Rx the mass 
and radius of the compact object and M its variable accre- 
tion rate. Assuming a neutron star, we will proceed adopting 
Mx = 1.4 Mq and Rx = 10 km. If the wind accretion onto 
the compact companion can be described as a spherical inflow 
(Bondi[1952), the accretion rate as a function of the binary sep- 
aration r is given by: 



M = 47t{GMx) 



^el 



(r) 



(2) 



where v„(r) is the wind velocity, p(r) is the density of the stellar 
wind and Vi-si(r) its supersonical velocity relative to the neutron 
star. We will approximate the wind velocity Vv,{r) by a typical 
jS-law appropriate for O-type stars: 

(3) 



Vw('-) 



1 - 



R, 



Here, Voo is the wind velocity at infinity and /3 a power index. 
Assuming a simple radial flow, the corresponding wind density 
can be derived from the continuity equation as: 



opt 



(4) 



Following McSwain & Gies <2002> . suitable values for 



LS 5039 are 



2850 km s"', /? = 0.8 and M, 



opt 



8 X 10" Mq yr" . To proceed, we have solved numerically the 
Keplerian motion of the neutron star around the Mopt - 40 M© 
companion calculating both the eccentric and true anomalies 
(E and V). The relative velocity between the wind and the neu- 
tron star is then given by: 

2 n2 

na 



V?el('-) 



■ sin/i + v„('')cos V 



■ cos E - v„{r) sin V 



(5) 



where a is the semimajor axis and n - 2n/Poib is the mean 
angular motion. 
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